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Abstract 

A new physical mechanism is suggested to explain the universal depletion of high meson exci- 
tations. It takes into account the appearance of holes inside the string world sheet due to qq pair 
creation when the length of the string exceeds the critical value Ri ~ 1.4 fm. It is argued that a 
delicate balance between large N^. loop suppression and a favorable gain in the action, produced 
by holes, creates a new metastable (predecay) stage with a renormalized string tension which 
now depends on the separation r. This results in smaller values of the slope of the radial Regge 
trajectories, in good agreement with the analysis of experimental data in Refs. Q. 
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I. INTRODUCTION 



Recently a number of radially excited mesons have been experimentally observed ]T[-@] 
and it was discovered that for excitations like aj{2P), U3{2D), and p{2D) their masses are 
100-200 MeV lower than the theoretical predictions in different models, in particular in the 
relativized potential model (RPM)[Q and in the flux tube model [^. A fC-matrix analysis 
of the Crystal Barrel data has shown that the Regge trajectories as a function of the radial 
quantum number Ur continue to be linear up to high excitations like the four L states 
= 0, 1, 2, 3) and can be described by the n^-trajectory (the radial Regge trajectory) [Q: 

M'^{nr,L) = M^{0,L) +nLnr (fixedL) (1.1) 

where the slope Ql was found to vary in the narrow range 1.15 < 0, < 1.30 GeV^ for different 
L-wave states. 

In Ref. [0] the orbital excitations of the light mesons (n^ = 0, L < 5) were studied 
in detail in the framework of the QCD string approach [|| and the spin- averaged meson 
masses M{nL), the Regge slope a^, and the intercept ^^(O) were calculated analytically 
and expressed through a single parameter, the string tension a, while the Regge intercept 
0^(0) does not depend on a and is a universal number. The calculated values of M{nL), 
a^, and ai^O) turn out to be in very good agreement with the experimental data. (Note 
our notation: a state denoted by nL has radial quantum number rir = n — 1, so the lowest 
state with a given angular momentum L is IL and has Ur = 0.) 

The situation appears to be different for the radial excitations (called also radials in what 
follows) calculated also for the linear string potential (or the linear plus Coulomb potential) 
with the same string tension. Thus the masses of the radials obtained do lie on a linear 
trajectory Eq. (|1 . 1| ) which has a slope that is practically independent of L, but its value 
VLq ^ 2.0 GeV^ appears to be a factor 1.6 - 1.5 larger than f2(exp), the value extracted 
from the experimental data (see Sect. 2). In particular, the masses of the second radial 
excitations (and even some first ones like 2D and 2F) are 100-150 MeV (n^ = 2) (50 - 100 
MeV for rir = 1) higher than M(exp) from Refs. - 0- 

This phenomenon, the lowering of the masses of the radials, was already discussed in 



Refs. and where it was supposed that this effect is connected with the opening 



of new channels, i.e. with hadronic shifts. However, hadronic shifts cannot produce the 
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global, quantum-number independent shift down of all radials on the n.f.-trajectory having 
been observed in experiment. In particular it cannot provide almost the same slope fl in 
Eq. ( |1 . 1| ) for different L, since hadronic shifts strongly depend on the quantum numbers of 
a decay channel, the closeness to the decay thresholds, the widths, and many other specific 
features of meson decays [^. 

Thus a basic paradox of mesonic spectra is that there are highly excited meson states 
with large width, implying strong coupling to decay channels, which nevertheless lie on 
linear Regge trajectories. This situation implies that first, the string between a quark and 
an antiquark continues to exist up to large energy excitations and can be as large as 2.5 — 2.8 
fm (see below). Secondly, quark pair creation does not dominate for such excitations, in 
particular the probability of string breaking is not large. How to reconcile these conclusions 
with strong decays (large width) of the Regge-string mesons? One can argue that pair 
creation is suppressed as l/N^ at large A^^^. Moreover, in experiment this parameter appears 
to be 1/10 rather than 1/3 for iVc = 3 which can be seen in the width to mass ratio for large 
excitations. The present paper suggests at least a partial answer to these questions taking 
as a characteristic example the radially excited mesons and generalizing to all highly excited 
light mesons which serve as a good illustration to the paradox stated above. 

So, to explain the "global correlated shift down" of the radial excitations we suggest 
here an alternative physical picture, which in first turn takes into account the behavior of 
the string in highly excited hadrons, and in addition the specific character of the vr-meson 
interaction with a light quark (antiquark) and the string connecting them. 

In contrast to hadronic decays like p — p, p — qq, u — u etc., which may occur due to qq 
loop creation inside the string's world sheet and subsequent string breaking, the 7r-meson, 
as well as other Goldstone particles, locally interacts with a quark (antiquark) sitting at 
the end of the string [|1^ and therefore the string may not break due to the emission of a 



TT-meson from a quark. 

The same statement is true when the creation of a vr-meson is accompanied by a p, ao 
or /o meson, since these mesons are actually described by the remaining string in the final 
state. Therefore it is natural to assume that below the threshold of p — p, {u — uj etc.), i.e. 

Ethr ^ 2Mp ±Tp> 1.4 GeV (1.2) 

the string effectively stays intact. When only channels like (n7r)7r, {mr)ri, it ~ p, r] — p etc. 
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are open, we argue that the string between a hght quark and antiquark is not modified by 
open channels and has the same string tension ctq ^ 0.18 GeV^ as the string between a static 
quark Q and antiquark Q. Then one can determine the characteristic size of the string Ri 
which corresponds to the value i^thr ~ 1.4 GeV: 

GqRi ~ Ethr,ori?i > 1.4GeV/ao = 1.45 fm. (1.3) 

Thus our first assumption here is that up to distances ~ 1.4 fm the string potential is not 
distorted by meson decays. This statement is in agreement with lattice calculations where 
in the presence of dynamical fermions the static potential appears to be the same as in the 
quenched approximation up to the separations of order 1.2 — 1.5 fm [|T^, | |14[ |. 



For higher excitations, E* > Ethr, the qq pair creation inside the world sheet (with the 
quantum numbers ^Si and ^Pq) is already possible and at first sight the problem is becoming 
essentially a two- or many-channel problem. However, we shall assume and argue here (see 



Section |T|) that up to very high excitations, Ecr > 2.5 GeV, i.e. in the range 

1.4 GeV <E*< E„ ~ 2.5 GeV (1.4) 

and for the time extensions 5T ~ l/F < 1.0 fm, only virtual loops or loops of small sizes 
are created. As a result the probability of hadronic decays like M^p — p, uj — uj, is small 
while with a large probability the string remains unbroken. This assumption is necessary 



to explain the linear character of the rir-trajectories Eq. (|1 . 1| ) up to high excitations of the 
order of 2.5 GeV. 

However, in the presence of such virtual or small qq loops the string tension is renormalized 
and becomes dependent on the separation r. In the suggested picture (when at small r, 
r < Ri, a = constant = o"o) the attenuation of the string tension is being felt only at 
distances r > Ri and continues at least up to the value, i?2 ~ 2.5 fm with o"(i?2) = <y2- It is 
important that this "asymptotic" value a2 strongly affects the slope fl of the n^-trajectory 
given in Eq. (|1 . 1|) . 

The most important feature of this picture is the existence of a prehadronization stage 
in the range of string parameters, 

Ri ~ 1.4 fm < r < i?2 ~ 2.5 fm, T <Tcr ^ 1.0 fm, (1.5) 

where the string tension depends on r while at the same time the string with large probability 
remains unbroken. The dynamics in this prehadronization region can be effectively described 



in the one-channel approximation taking into account virtual quark loops and the open 
hadronic channels (mostly like {mi)7i, n — p, rj — p etc.) through a universal dependence of 
the string tension on r: 

a = constant = cxo, r < Ri, 

a = a(r) Ri < r < R2 (1.6) 

At the present stage of the theory the function a{r) is not yet calculated in full QCD and 
therefore we formulate here the problem in a different way: how to extract information 
about o"(r), or the string breaking, from highly excited meson masses, in particular from the 
slope of the radial Regge trajectories. 

We shall show here that there exists a direct connection between the slope Q and the 
two most important features of a{r): the value of Ri where the string tension is becoming 
r-dependent, and the value (T2 which characterizes the string tension in the region where 
breaking is already possible. In lattice calculations a flattening of the static potential due to 



qq pair creation at r ~ 1.2 — 1.5 fm seems to be observed |jTj], unfortunately, lattice points 
have very large errors and at present definite conclusions about the exact value and form of 
the static potential at large r cannot be derived from lattice measurements. 

We concentrate below on these considerations and suggest a workable and simple model 
for the mesons of large size, both radial and orbital ones, which yields meson masses in 
good agreement with experiment. It will be shown that in the proposed picture the linear 
n,f.-trajectories with a rather small slope, Q ~ 1.3 GeV^ (~ 1.5 GeV^ for the spin averaged 
S-wave states) close to the experimental numbers, can be easily obtained. 

The plan of the paper is as follows. In Section ^ the analytic formulas for the Regge 
slope and the intercept for the linear potential are derived and the masses will be expressed 
through a single scale parameter - the string tension. It will be shown that for the standard 



linear potential the slope fig is a factor 1.6 larger than in experiment. In Section |T| the 
effects of gg-pair creation (unquenched situation) on the meson masses of large radii are 
discussed and in Section |V| a modified nonperturbative potential is proposed for which 
the meson masses are calculated. In Section the Regge slope and intercept of the n^- 



trajectories are presented. In section our Conclusions and some prospectives are briefly 
discussed. In the Appendix the results of the detailed calculations for the meson spectra are 
included. 



5 



II. THE PROBLEM OF HIGH EXCITATIONS FOR THE LINEAR POTENTIAL 



In the QCD string approach the Hamiltonian is derived from QCD under definite and 
verifyable assumptions. First, the times and distances involved are considered to be larger 
than the gluonic correlation length Tg-. r ^ Tg {Tg ^ 0.2 fm in lattice calculations ||15|). This 
condition is always valid for the light mesons having large sizes i? > 0.8 fm. Secondly, the 
string (hybrid) excitation scale is large, AMgtr > 1 GeV, and therefore in first approximation 
the meson and hybrid excitations are disconnected. 

Then the spin-averaged mass M{nL) of the light meson with arbitrary quantum numbers 
nL {L < 5)are determined by the following mass formula 0, ^: 

M{nL) = Mo{nL) + Antrim) + AsE{nL) (2.1) 

where Mq is the eigenvalue of the unperturbed string Hamiltonian H}^ coinciding with the 
spinless Salpeter equation (SSE): 

2y5^Tm2 + V{r)) ip{nL) = Mo{nL)i){nL), (2.2) 



where m is the current quark mass taken here to be equal to zero (m = 0). The potential 
V{r) contains in the general case both a perturbative part, the Coulomb interaction Vcir), 
and the nonperturbative string potential ^p(r) with a string tension that is in general 
dependent on the separation r 0: 

An 

V{r) = --^ + a{r)r. (2.3) 

It is instructive to consider first a linear potential with a = constant = ctq with a mass 
formula that is more transparent and can be presented in analytical form. In this case the 
string correction was defined in [0, |l^: 

As.(nL) = _-o<r-\>L{L + l) ^ _2a, < r-^ > LjL + I) 

^ ' 8fil{nL) M^inL) ' ^ ^ 

where the following relations valid for the potential o"or are used: 



MoinL) = 4fio{nL), (aor) = 2/io(nL). (2.5) 

The constituent mass fio{nL) was derived to be the average of the quark kinetic energy 
operator 0,11: 

/io(nL) = {.JfP + m^)nL. (2.6) 



From the definition Eq. ( p.6|) it is clear that /iq depends on the quantum numbers nL of a 
given state and can be expected to grow for high excitations. 

An important contribution to the meson mass Eq. ( |2.1| ) comes from the nonperturbative 
self-energy term Ase('^-^^) calculated in Ref. |T^: 

AsE(nLj {!.() 

TT/io TTMQ[nL) 



where the parameter ri{f) depends on the quark flavor / and can be calculated, see Ref. [jlT 



ri{nn) = 0.90 (2.8) 

Here it is worth to notice that the Coulomb corrections Ec{nL) to the light meson masses 
are small, \Ec\ < 100 MeV and can be neglected in the first approximation. The results 
of the exact calculations for a linear plus Coulomb potential Eq. (|2.3| ) will be presented in 



Section |^ and in the Appendix. 

The self-energy term enters the squared mass in such a way that it gives rise to an 



important negative constant Cq, 



Co = (2.9) 



TT 

SO that the mass formula for the squared meson mass is 

i.e. it is proportional to a and does not contain any free parameter. In Eq. (|2.10|) the con- 
tributions of the small terms A^^^. and 2AstrAsE were neglected to show the most important 
features of the meson spectra. On the contrary the term Ag^^ > 0.10 GeV^ is not small for 
any state and therefore is kept in the mass formula (|2.10| ). 



An important next step refers to the approximation for the eigenvalues M^{nL) of the 
SSE valid for the linear potential dor : 

Mq (approx) = [8L + 47r^(nL) rir + 37i]ao. (2.11) 

This formula with ^ = 1.0 reproduces the exact values of Mq(?2L) with an accuracy better 
than 2% for all riL-states with the exception of the IS* and IP states where the accuracy 
is 3-6% [0. In Eq. ( p.ll|) the coefficient ^{nL) ~ 1.0 weakly depends on L and slightly 



decreases with growing n^, e.g. for the 45" state ^(45) = 0.99 while ^(4F) = 0.96. In what 
follows in most cases we put i{nL) = 1.0. 



Then with the use of the expression ( |2.11| ) and redefining the matrix element < >= 
y/oQ < > where < p^^ > is aheady independent of ctq, one can rewrite the mass formula 
Eq. ( p.lO| ) in the following way, 

M^(nL) = {a'^y^ L + [Anrir + b{nL)](To (2.12) 

giving for the L-trajectory, 

L = a'^ M^{nL) + ^^(n), (2.13) 
the Regge slope (in general 7^ 0): 

^2< p-^ > {L + l) 



L 



2 8 



(2.14) 



For the leading L-trajectory (< p ^ > a/L + 1 ~ constant) the slope is given by the constant 
i 



or 



1 



a? 



6.95ao 



(6.95±0.02)ao 
O.SOGeV-^ for (To = 0.18GeV^ 



(2.15) 



which is in good agreement with the experimental value a^(exp) = 0.81 ± 0.02 GeV ^. For 
the orbital excitations with a fixed Ur > 0, a'i{nr) appears to be a bit smaller since for them 
6{nL) in Eq. {^J§ is smaller, e.g. 5(n^ = 1) = 0.85 ± 0.05 and (5(n^ = 3) ~ 0.50 while 
S(nr = 0) = 1.05. Then for ctq = 0.18 GeV^ 

a^(n, = 1) = l/[7.2(To] = 0.77GeV-^ 

a'^(nr = 3) = l/[7.5(To] = 0.74GeV-^ (2.16) 

However, this slope gives larger values for the radial excitations, e.g. M{2P) = 1.82 ± 0.03 

GeV while the expected experimental number is M{2P) < 1.70 — 1.75 GeV 0, 0. 

From Eq. ( p.lQ ) it follows that the constant b{nL) in Eq. (|2.12| ) is 

32r] ?,2rf 



h{nL) = 37r 



(2.17) 



TT 7r2(L + ^ + f)" 
As was stressed in Ref. 0], for the h{lS) meson it is more precise to use the exact eigenvalue 
with Mq{1S) = 9.967 ctq instead of the approximation ( p.llD so that for the leading L- 
trajectory [i] = 0.90) one finds 



TT 



0.257 + 



2567^2 
9.967 7r2 



2.365. 



(2.18) 



Then from the mass formula ( p.l2| ) the intercept of the leading L-trajectory, 

aLiM^ = 0) = -a'^bilS) do = -2.365/6.95 = -0.34, (2.19) 

does not depend on the string tension and is a universal number. With 10% accuracy it 
coincides with the experimental number aL(0)exp = 0.30 ± 0.02 0. 
For rir the intercept of the L-trajectory is 

aL{nr,M'^ = 0) = — a'^(nr) [4?™^ + 6(nL)]cro 

0.257 + ^fJi^ ,^ + ATTUr 

= S-f{nL) ^""^ ^ ^^-^^^ 

In the intercept ( |2.20| ) the term b{nL) is small as compared to Ann (the largest correction is 
for the S'-wave states but even then it is < 10%) and can be neglected in first approximation. 
Then also neglecting 6{nL) compared to 8 one obtains 

aLirir) - (2.21) 

which does not depend on L and cxo. For the neighbouring intercepts the difference is equal 
to the following constant, 

aL{nr + 1) - aiirir) = -n/2 = -1.57. (2.22) 

The magnitude of the intercept is rather large due to the presence of the large number 
47rnrO"o in the eigenvalue MqIuL) (or ( 2.12 )). Note that another large number, Snao, 



which is present in MqIhL), is practically cancelled by the constant Cq, Eq. ( p.9| ), coming 
from the self-energy contribution. Now one can present the masses of the radials in the form 
of the n^-trajectory Eq. ( |1 . 1| ) : 

M'^{nr, L) = {0, L) + no rir (L fixed), (2.23) 

where from the mass formula (|2.10| ) one finds for the linear potential the following expression 
for 

no = [47rao^(nL)][l - xinL)] ~ 1.90 - 2.10 GeV^ forao = 0.18 GeVl (2.24) 

(xinL) = ibinL) ~b{n + l L))/Att) 

In Eq. (|]24D 0.95 < ^{nL) < 1.05 and 0.05 < x < 0.10 for all nL states {ur < 4, L < 4) 
so that VLq is only 5%-10% smaller than AnaQ. 



Thus for the hnear potential we have obtained 

(i) The masses of the radials he on hnear rir-trajectories as in Eq. ( |1 . 1| ) ; 

(ii) The slope Qq does not depend on L (with 95% accuracy) in agreement with experimental 
observations; 

(iii) However, numerically the value of Qq turns out to be ~ 1.6 times larger than the value 
extracted from the experimental data ||^. 

To explain this phenomenon we suggest below a physical mechanism which can be applied 
to highly excited mesons. 



III. MESON MASSES AND QUARK PAIR CREATION 

The effective Hamiltonian derived from the QCD Lagrangian with the use of the Fock- 
Feynman-Schwinger representation [p!^] is based on the quenched approximation, where the 
quark determinant is replaced by unity. Based on the width-to-mass ratio and on the 
existence of linear Regge trajectories for the mesons it is usually argued that the effects of 
the sea quark loops coming from the quark determinant cannot be large and are estimated 
to be around 10%. The same estimate of this correction is obtained from lattice calculations 



for the unquenched low-lying hadrons[|l^]. Moreover, lattice calculations of the QQ static 
potential up to 1.0 — 1.5 fm do not show a significant difference between quenched and 
unquenched calculations ||T3[|- [p0[ . 



It can be shown that the radial excitations and high orbital excitations {L > 4) have 
sizes exceeding 1.5 fm (see Table |) and therefore one should reconsider possible effects 
of quark loops on the large-size mesons. A dedicated study on the lattice ]14[ shows a 



flattening of the static potential or decreasing of the string tension at separations r > 1 fm. 
However, at present lattice points have large errors at such distances, quickly deteriorating 
with increasing r and one cannot extract the exact form of the static potential at large r 
from lattice data. 

Above the threshold of qq pair creation the static QQ pair could decay into two heavy- 
light mesons with mass Mhl = 2mQ + 2i?HL where i?HL is the excitation energy of the 
heavy-light meson which can be calculated in the framework of the formalism presented in 
Ref. and is found to be EuLi^s = 0) = 0.73 GeV and Euhias = 0.39) = 0.53 GeV for a 



6-quark mass rrib = 4.8 GeV. 
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For the c quark with rric = 1.4 GeV the values are close: i?HL(«s = 0) = 0.76 GeV and 



Euhicts = 0.39) = 0.58 GeV are obtained|^. So one can expect that for lighter quarks i?LL 
is also ^ 0.7 — 0.6 GeV and the threshold of qq pair creation inside the world-sheet of the 
string to be Mthr = 1-2 — 1.4 GeV in accord with the experimental p — p threshold Eq. ( |1.2[ ). 
Expressing this value in terms of the distance Ri: Mthr = <^oRi one finds the number 

Ri = 1.3 - 1.5 fm 



close to the value Ri in Eq. ( |1.3|) . 

At this point one should stress that the phenomenon discussed, viz the pair creation just 
on the string, does not necessarily exhausts all possible meson decay mechanisms or the true 
hadronization of the mesons. Namely, as was shown in Ref. |T2[, pions are directly coupled 
to the quark (antiquark) at the ends of the string and can be emitted from there without 
breaking it. 



Indeed, bosonization of quark degrees of freedom in Ref. leads to the following term 
in the Lagrangian in the local limit: 



A/:^^) = J dtSx 



g(x)a|f|75^^g(x) 



(3.1) 



where the string starts from the antiquark position x = 0; the field vr" is the octet of 
Nambu-Goldstone mesons, and F^^ = 93 MeV. 

The same operator A^*^^-' between quark bound states can be rewritten with the use of 
the Dirac equation as follows 

A/:« = (7^Trg7'^75Cu^g, (3.2) 

with 

^ti = (ud^u^ - u^d^u) u = exp (^^75^^^ , 9a = l-O- (3.3) 

From this expression one can see that AC^^^ describes the emission of an arbitrary number 
of Nambu-Goldstone mesons from the quark position and therefore may describe pionic and 
double pionic hadron decays while the string plays the role of a spectator and stays intact. 

At the same time qq pairs around the string (sea quarks) should be identified with the 
loops of the determinant which can be written as ||2^, |23[ : 

d^xj^ -j-exp{-m^s){Dz),,expi-J^^\w{C). (3.4) 
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The integral {Dz)xx in Eq- ( |3.4| ) is taken along the closed loop C from some point x back to 
X and contains an integral over loops of all sizes. Alternatively, one can separate the quark 
determinant into parts of small and large eigenvalues as it Ref. [Q, 

det{A) = det ir(v4) det uv(^), (3.5) 

where detiR takes into account the small eigenvalues (large loops) A„ < Acut while detuv 
contains the large eigenvalues. The latter correspond to the contribution of small virtual 
loops which should be properly renormalized. 

In the physical picture detiR corresponds to the chiral effects which disappear in det(m^ + 
D"^) for large m^: ^ ^qcd together with effects of large loops. 

At this point it is important to stress that treating the resonances as quasistationary 
states one should always consider Green's functions and correspondingly Wilson qq loops of 
finite time extension T, T ~ l/r(nL) ~ 1 — 2 fm. 

For such finite times one can write a general expansion for the original Wilson loop of 
the qq meson with pair creation due to the quark determinant Eq. ( |3.4] ) as follows 



{W{C) det(m + D)) = Wo{C) + ^Wi{C, Ci) + ^W^2(C, Ci, C2) + . . . , (3.6) 

where the coefficients Oj = 0(1) and higher Wilson loops Wi,i > 1, are averaged products 
of the i + 1 Wilson loops, i.e. 

W,{C,C^) = {{W{C)WiC,))), WiC) = ^Tr#(C). (3.7) 

It is clear that when T — * 00 and r > Ri (so that decay is energetically possible) the 
asymptotics of the r.h.s. of Eq. (|3.6|) is 



exp(-arT) + exp(-T(Mi + M2)) + ^ exp(-T^ M,) + . . . , (3.8) 

where J2 is the sum of the masses of the decay products {i = 1, 2 in the simplest case and 
i = 1 . . .k for meson decay into fc-particles etc.). So, if T — > 00 the resonances are dying 
out and the second term in Eq. (|3.8| ) is dominant, which means that only the products of 
meson decay are left. 

However, for finite T ^ 1 — 2 fm the situation is different and one can expect a delicate 
balance between the large Nc limit (suppressing qq pair creation) and large times T (pre- 
ferring large internal loops). This statement is the dynamical basis of our main assumption 
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about the existence of a specific state of the string with effectively a hole inside the world 
sheet which obeys the area law with a reduced (renormalized) string tension a* = cr(r). 

Thus we introduce the new concept of a transitional regime which in terms of QQ sepa- 
rations r refers to the region 

Ri<r <R2 

where Ri ~ 1.2 — lA fm. R2 ~ 2.5 fm corresponds to high excitations with energies E* ~ 2.5 
GeV. 

In this region, due to relatively small quark loops and correspondingly small holes in the 
string world sheet, the string tension is nonperturbatively renormalized and decreases with 
growing r. While the loops are still virtual, their presence does not lead to actual string 
breaking (with appreciable probability) even though the energy of the string for r > i?i is 
sufficient for meson decay. We suggest here that only for much larger distances, r > R2 — 2.5 
fm the string breaking happens with large probability. 

A first argument in favour of this picture is that highly excited mesons of large radii 
do exist and lie on the corresponding linear Regge trajectories, while their characteristics 
can be computed as in the QCD string approach, neglecting decay channels. At the same 
time at such large distances the string tension cannot remain intact and should be strongly 
decreased by the appearance of holes inside the string world sheet due to pair creation. 

The second argument refers to the experimental information about strong decay modes. 
The first strong decays without Goldstone particles is observed only for /o(1370) ^ p — p 
and /2(1565) — > po ~ Po decays in accord with our picture that a = constant at distances 
r < Ri or excitations < 1.4 GeV. 

For higher excitations, in the range 1.4 GeV < E* < 2.5 GeV, at present only several 
strong decays like /2(1640) — > c<j — c<j, (seen); 7r2(1670) uj — p (branching ~ 2.7%), and 
/4(2050) uj—uj (branching ~ 25%) have been measured and in all these cases the branching 
ratios of decays without Goldstone particles are never large. 

The third argument refers to the r.m.s. radii R{nL) = V< >nL of the radials which 
are (n^ 7^ 0) calculated for the linear potential with ctq = 0.19 GeV^ (see Table |). (It is 
worth to notice that the numbers given for R{nL) represent the lower limits of the true 
r.m.s. radii, since they correspond to larger meson masses and the actual values of R{nL) 
are about 20-50% larger (see Table [V^ ). From the values of R{nL) given here one can see 
that among the ground states only the IG state has R{1G) ~ 1.4 fm, while for the 2L and 
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TABLE I: The light meson r.m.s. radii R{nL) in fm for the a^r potential with do = 0.19 GeV^. 



nj.\L 





1 


2 


3 


4 





0.80 


1.02 


1.21 


1.37 


1.51 


1 


1.22 


1.31 


1.52 


1.70 


1.77 


2 


1.53 


1.66 


1.78 


1.90 


2.00 



the 3L states (with the exception of the 25* state with relatively small r.m.s. radius equal 
to 1.21 fm) i?(nL)are in the range: 



1.4 fm < R{2L) < 1.8 fm (n^ = 1), 
1.7fm < R{3L) < 2.0 fm (n^ = 2). 



(3.9) 



The values of M(exp) for all states underlined in Table | are shifted down compared to the 
theoretical values calculated with the same linear potential that gives a good description of 
the orbital excitations with 71^ = 0. This example agrees with our estimate of the charac- 
teristic size Ri = 1.4 fm where the pair creation is beginning to affect the string tension 
a(r). 

The r.m.s. radii in Table | also show that even for the linear potential R{nL) > 1.90 fm 
for such states as the 3F states: ^(2290), /3(2280), a4(2280), and a3(2310)||. 

At this point it is important to stress the difference and similarity of our approach with 
that in Refs. |1^. In both approaches it is stressed that the qq pair creation is responsible 
for the renormalization of the string tension and therefore the unquenched string tension is 



lower than in the quenched case. Moreover in Ref. |T0| as well as in the present paper it is 
emphasized that even when level crossing occurs, i.e. when V{r) equals 2Mhl, the string 
potential can be used throughout, only the string tension being renormalized. 

The difference between both approaches is in the meaning of this renormalization. We 
assume here that there exists a universal (quantum number independent) prehadronization 
stage when small quark loops attenuate the string tension, while one can still neglect the 
influence of specific decay channels which produce the hadronic shift. 

It was already realized in Refs. |10| that the very fact of the occurrence of smooth 
Regge trajectories and ordered hadronic spectra is difficult to explain if the hadronic shifts 
are essentially important, since the latter depend on the concrete hadronic channels involved 
and vary irregularly from channel to channel. 
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In the suggested picture due to the universal predecay (prehadronization) stage with 
renormahzed (attenuated ) string tension not only can the spectrum be calculated, but the 
notion of linear Regge trajectories, both radial and orbital, is kept intact. 



IV. MODIFIED CONFINING POTENTIAL 



From the physical picture discussed in Section |ITl| and the introduction it follows that up 
to the characteristic distance Ri ~ 1.2 — 1.4 fm the string tension a is constant, while for 
larger r it depends on the qq separation as in Eq. (|1.6|). We propose here the nonperturbative 



(4.1) 



potential a{r) r with the string tension taken in the following form: 



Vnp (r) = a(r)r, a(r) = (Tq 



1-7- 



exp(v/o'o(r - Ri)) 



5 + exp(v/^Q(r-i?i))_ 

In the definition Eq. ( [4.1| ) the constant 7 determines the value of the string tension at large 
separations so that 

= (T(r > 2.5 fm) ~ (To(l -7). (4.2) 

Note that do ~ 0.18 — 0.19 GeV^ defines the common scale of the modified string potential 
and can be fixed by the Regge slope of the leading L-trajectory. 

In general our calculations are performed with the potential V{r) Eq. ( |2.3| ) containing 
the perturbative potential Vp{r), 

VAr) = (4.3) 

where for the strong coupling constant the value a.^ = 0.30 is taken. 

The best description of the meson spectra was obtained for the following set of the 
parameters in a{r): 

7 = 0.40, Ri = 6 GeV-\ B = 20.0, (4.4) 

while the values of (Jq = 0.185 ±0.005 GeV^ and = 0.30 ±0.08 can vary in narrow ranges. 
In Fig. this potential is drawn for the parameters Eq. ([4.4| ), = and ctq = 0.19 GeV^. 

At the distance R2{nL) = 2.5 fm the value (J2 = 0.116 GeV^ turns out to be rather small, 
a value 40% smaller than ctq = 0.19 GeV^. 

The spin-averaged meson masses M{nL) are calculated solving the SSE Eq. (|2.2| ) with 
the modified potential V{r). Their values will be given in the next section. Here only some 
characteristic features of the SSE solutions for the modified potential will be discussed. 
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r (GeV"') 

FIG. 1: Modified potential with parameters given in Eq. ( [4.4[ ). For reference a simple linear 
potential with o" = 0.19 is also plotted. 

TABLE 11: The values of a{nL) (in GeV^) for the potential V{r) Eq. ( |4.1| ) with the parameters 
Eq. (U) (o-o = 0.19 GeV^, = 0.30). 



nr\ L 





1 


2 


3 


4 





0.188 


0.187 


0.184 


0.179 


0.168 


1 


0.183 


0.176 


0.164 


0.153 


0.150 


2 


0.165 


0.155 


0.154 


0.154 


0.150 


3 


0.156 


0.156 


0.154 


0.154 


0.150 



First, the average values of the string tension a =< (j{r) >nL turns out to be almost 
constant for the ground states (n^ = 0, L < 3), while for the radials with > 2 d" is already 
20% smaller, see Table ||. 

Secondly, the kinetic energy and therefore the constituent quark mass fi{nL) does prac- 
tically not change for the ground states in contrast to the high excitations where the string 
size is growing and the quark kinetic energy is becoming significantly smaller. In Table |T| 
fi{nP) for the nP states for the potential Eq. ( [4.1| ) with = and the linear potential a^r 
are compared (in both cases do = 0.182 GeV^). 
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TABLE III: The constituent masses ^{nP) for the modified potential Eq. (4.1) with 
fj.o{nP) for the potential aor with (Jq = 0.182 GeV^ in both cases. 



and 








1 


2 


3 


4 


li{nP) 


0.424 


0.430 


0.441 


0.559 


0.616 


lioinP) 


0.451 


0.582 


0.697 


0.787 


0.872 



TABLE IV: The spin-averaged meson masses of the nS states (in GeV) for the modified potential 

0.30 and gq = 0.19 GeV^ *. 



V{r) Eq. (O) with the parameters Eq. 



nS 


IS 


2S 


3S 


45 


55 


M{nS) 


0.618 


1.400 


1.868 


2.176 


2.502 




(0.673) 


(1.520) 


(2.122) 


(2.602) 


(3.006) 




0.612 


1.41 ±0.02 


7r(1.80) 


p(2.15) 





The numbers in brackets correspond to the Corne 



1 potential Eq. (|5.1|) with the same ag = 0.30 



and do = 0.19 GeV^ 



From Table |TTl| one can see that for the 4P(5P) states the difference between the con- 
stituent masses is large and reaches ~ 30% . Note that for the modified potential the 
relations (|2.5|) are not valid anymore and 



(4.5) 



and therefore in this case the string correction Agtj. is given by a more general expression 
than in Eq. (U) 0: 

= _i4llii(£±i). (4.6) 



/i(6/x+ < a{r)r >) 

Also in the self-energy term AsE{nL) Eq. ( p.7|) , a must be used instead of ctq- 



V. THE MASSES OF THE RADIALS 

The masses of the radial excitations for the modified potential Eq. ( ^.1| ) with the Coulomb 
interaction included are presented in Tables P^[VI| for the parameters given in Eq. ( [4. 4] ), 
tts = 0.30, and the string tension ctq = 0.19 GeV^. In all cases the calculated meson masses 
turn out to be in good agreement with the existing experimental data. 
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TABLE V: The spin-averaged meson masses of the nP-states (in GeV) for the potential V(r) 
Eq. (O) with the parameters Eq. (O), = 0.30 and uq = 0.19 GeV^ * 



nP 


IP 


2P 


3P 


4P 


M{nP) 


1.190 


1.715 


2.090 


2.388 




(1.263) 


(1.933) 


(2.438) 


(2.859) 




1.252 forM(a,/(lP)) 


ai(1.70) 


ai(2.10) 


ai(2.34) 






02(1.75) 


ao(2.05) 


/o(2.34) 




1.245 for M(/j(lP)) 


/2(1.65) 


/o(2.095) 




*) see footnote to Table 


V. 



The numbers in brackets in Table are the masses calculated for the Cornell potential, 

4a, 



Vrir) 



3r 



+ (Tor. 



(5.1) 



From a comparison of the numbers given one can see that with the modified nonperturbative 
potential the mass of the 43(53) states appears to be 400 MeV (500 MeV) lower than for 
the Cornell potential Eq. (|5.1|) , while for the IS* and IP states the difference is only about 
60 MeV. 

We observe the same picture for the P-wave and higher excitations when already for 
the 2P (3P) states the spin-averaged mass is ~ 200MeV (350 MeV) smaller due to the 
modification of the string potential (see Table |V|, where the numbers in the parentheses are 
calculated with the Cornell potential Eq. ( |5.1| ) with the same a, and ctq). 

The calculated masses of the n3 and nP states as well as the nD and nF states (see 
Table |VJ) appear to be in good agreement with experiment. However, for the nD and nF 
states a better agreement is obtained for smaller values of the strong coupling constant and 
the numbers given in Table ^ refer to a, = 0.21. This fact may be connected with a 



suppression of one-gluon-exchange for large-size mesons. 

Thus one can conclude that due to the attenuation of the string tension in the potential 
Eq. ( |4.1|) the masses of the radials turn out to be ~ 100 — 200 MeV (for the 2L states), 
~ 300 MeV (for the 3L states), and ~ 350 — 400 MeV smaller than for the standard linear 
potential. It is of interest also to compare the r.m.s. radii R{nL) of the radials for the 
modified potential (see Table [VI1| ) with that for the aor potential given in Table |. 

From this comparison one can find out that the size of the 2P mesons is changing from 
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TABLE VI: The spin-averaged meson masses (in GeV) for the nD and nF states for the potential 
V{r) Eq. (13) with the parameters Eq. (O), = 0.21, and ctq = 0.19 GeV^. 



n 


ID 


2D 


31) 


IF 


2F 


3F 


M{nD) 


1.628 


1.973 


2.290 


1.926 


2.214 


2.480 


M(exp) 


P^2(1.67) 


vr2(2.0) 


7r(2.25) 


04(2.01) 


04(2.26) 






P3(1.69) 


P3(1.98) 


P3(2.30) 


03(2.03) 


03(2.28) 





TABLE VII: The r.m.s. radii R{nL) (in fm) of the nL states for the potential V{r) (41) with the 
parameters (43), = 0.30, and uq = 0.19 GeV^. 



state \L 





1 


2 


3 


4 


IL 


0.74 


1.03 


1.29 


1.61 


2.10 


2L 


1.32 


1.74 


2.22 


2.61 


2.74 


3L 


2.22 


2.53 


2.63 


2.67 


2.81 


4L 


2.58 


2.66 


2.81 


2.97 


3.09 



R{2P) = 1.31 fm to 1.74 fm while the mass M(2P) is shifted down by about 200 MeV. 



From Table |VII| it is also seen that for the modified potential highly excited radials, like the 
3P mesons, have very large r.m.s. radii ^ 2.5 — 2.8 fm, in particular the experimentally 
observed mesons p(4S') and aj{3P) have R{nL) ~ 2.5 — 2.6 fm. 



VI. THE SLOPE OF THE n^-TRAJECTORIES 



There are not many radial excitations with well established masses which are included in 
the PDG compilation ^ . Most radials were observed in the BNL and Crystal Barrel exper- 
iments and discussed in many papers for the last five years (see Ref. 0-0 and references 
therein). Here we present the values of the slope VL defining the n^-trajectory Eq. (|1 . 1| ) . 

Since we have calculated here only the spin-averaged masses of the radials, correspond- 
ingly just for them the n^-trajectory Eq. (|1 . 1|) will be calculated below. Although in many 



cases there exists a large uncertainty in the values of M{nL) we give below in Table VIII 
several well established masses taking into account that the spin splittings are small. 

Then taking the difference between the neighbouring values one can calculate the 
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TABLE VIII: The experimental spin-averaged masses of the radials 



L 


IS 


25 


IP 


2P 


3P 


4P 


M{nL) 


0.612 


1.42±0.04 


1.25±0.05 


1.70±0.05 


2.07±0.03 


~ 2.34 


L 


ID 


2D 


3L> 


IF 


2F 




M{nL) 


1.67±0.02 


2.00±0.02 


~ 2.30 


2.02±0.01 


2.30±0.01 





values of VL in Eq. (pTTI) : 



^exp(*S') 
^cxp 



1.64 ± 0.11 GeV^ 
1.35 ± 0.25 GeV^ 
1.21 ±0.16 GeV^ 
1.21 ±0.08 GeV^ 



(6.1) 



which in some cases have a rather large experimental error, but for L 7^ practically coincide 
with the value Vl = 1.15 — 1.30 GeV^ obtained in Ref. 0. Unfortunately, among the nS 
states the very important ones p{3S) and 7r(4S') are still not observed and the known value 
of Qcxp{nS) obtained from the difference of M^(2S') and M^(IS') appears to be 20% larger 
than in Ref. P| and close to our number Eq. (|6.3|) . 

From the meson masses given in Tables |rV| - |V| one can calculate now the theoretical values 
of the slope Qthr taking for the mass M^(IL) with ra^ = in Eq. ( |1 . 1|) the experimental 
number. Then one finds 

rithiP) = 1.38 ± 0.05 GeV^ 
QthP) = 1.29 ± 0.06 GeV^ 

Vltk{F) = 1.22±0.13GeV^ (6.2) 

However, for the nS radials the slope was found to be slightly dependent on rir- 

nti,{S) = 1.60 - 1.45 forn^ = 0, 1, 2, (6.3) 

which is close to the experimental value f2exp = 1-6 ±0.1 GeV^ obtained for the first excited 
state, while for higher S'-excitations Qth was found to be ~ 15% smaller and equal to 



M^{4:S) - M'^{3S) = 2.18^ - 1.87^ = 1.26 GeV^ 



(6.4) 



20 



We can conclude that in the physical picture where the confining potential is modified due 
to qq pair creation, the slope Vt is decreasing from a value ~ 2.0 GeV^ for the standard linear 
potential to values in the range 1.2 — 1.35 GeV^ for mesons with L 7^ 0. 

VII. CONCLUSIONS 

We have considered here the light meson orbital and radial excitations using the effective 
Hamiltonian derived from QCD under definite and verifyable assumptions. In the QCD 
string approach the spin-averaged meson mass M{nL) can be calculated through the only 
scale parameter - the string tension and does not contain any arbitrary subtraction constants, 



since nonperturbative quark mass renormalization is taken into account as in Ref . [0 . The 
suggested formalism allows to resolve three old painstaking problems: 

To determine the origin of the constituent mass for a light quark which is derived to be 
the average of the quark kinetic energy operator and can be computed through the string 
tension; 

To obtain the correct slope of the Regge trajectory when the string moment of inertia is 
taken in account; 

To obtain the correct absolute values of the light meson masses and as a consequence 
the correct value of the L-trajectory intercept (which refers to the spin averaged meson 
masses). The use of the L-trajectories is very convenient since they are universal, i.e. in the 
closed- channel approximation they are the same for isovector and isoscalar mesons. 

In Ref. this formalism was successfully applied to the orbital excitations with rij. = 
when for the linear confining potential the string tension was taken constant. However, in 
an attempt to describe the radial excitations one encounters a serious problem - the Regge 
slope of the n^-trajectories calculated with the same potential appears to be 1.5-1.7 times 
larger than in experiment. 

This phenomenon, the lowering of the masses of the highly excited mesons, is connected 
in our picture with large sizes of the high excitations, which can be as large as 2.5 fm and 
lead to the formulation of the concept of the predecay (prehadronization) region where due 
to qq pair creation the string tension is attenuated at separations r > Ri = 1.2 fm — 1.4 fm. 
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In this physical picture it is important to take into account the specific character of the vr- 
meson interaction with a hght quark which occurs at the end of the string and therefore the 
string does not break due to the vr-meson emission and this fact reconciles a high probability 
of pionic exchanges for the mesons with the existence of linear Regge trajectories. 

The explicit and very simple model of the modified confining potential where the string 
tension depends on the separation r for r > Ri, allows to obtain the masses of the radi- 
als in good agreement with experiment and may be considered as an explanation for the 
observation that the masses of high excitations are lowered. In particular the centers of 
gravity of the 2Pj and 3Pj multiplets appear to be lower by ^ 200 MeV and ~ 350 MeV 
respectively. The slope of the nj,-trajectory Q {L ^ 0) is found to be 1.35- 1.22 GeV^ in 
agreement with the analysis in Ref . |0 . For the nS states the calculated slope is found to be 
larger, ~ 1.5 GeV^. The mechanism of reduced string tension has a universal character 
and does not depend on the quantum numbers and concrete positions of open thresholds in 
meson decays. 

The authors would like to express their gratitude to the theory group of Thomas Jefferson 
National Accelerator Facility (TJNAF)for their hospitality. One of the authors (Yu.A.S.) 
was supported by DOE contract DE-AC05-84ER40150 under which SURA operates the 
TJNAF. This work was partly supported by the RFFI grant 00-02-17836 and INTAS grant 
00-00110. 

APPENDIX A: MATRIX ELEMENTS FOR THE SOLUTIONS USING THE 
MODIFIED LINEAR POTENTIAL 

Here we present some characteristics of the SSE solutions Eq. (|2.2| ) for the modified confin- 
ing potential Eq. ( |4.1| ) needed to calculate the spin-averaged meson masses. The eigenvalues 
of the Eq. ( |2.2| ) for the potential Eq. (^Tl]) plus Coulomb potential with parameters Eq. ( |4.4| ) 
and (To = 0.19 GeV^ are given in Table |^ for the S- and P-wave mesons (a^ = 0.30), and 
in Table ^ for the D- and F-mesons (a<j = 0.21). 

We give also the constituent masses and the matrix elements < r^^ > entering the string 
corrections A^trinL) and AssinL) while the average values of a =< cr(r) >„,l are given in 
Table 0. 
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TABLE IX: The eigenvalues of the SSE Eq. (|2.2| ) in GeV for the Coulomb plus modified confining 
potential Eq. (O) with the parameters Eq. (O), ctq = 0.19 GeV^ and = 0.30 (L < 4, n,. < 3) 



nr\ L 





1 


2 


3 


4 





1.204 


1 721 


2.098 


2.402 


2.645 


1 


1.858 


2.190 


2.436 


2.628 


2.799 


2 


2.268 


2.471 


2.656 


2.840 


3.013 


3 


2.542 


2.730 


2.904 


3.065 


3.218 



TABLE X: The constituent masses (in GeV) and the matrix elements < r ^ > (in GeV ^) for the 



P- and D-waves for the same potential as in Table IX 





IP 


2P 


3P 


4P 


< r^^ > 


0.464 
0.250 


0.484 
0.203 


0.535 
0.180 


0.587 
0.168 




ID 


2D 


3D 


4L> 


fi{nD) 

< r^^ > 


0.526 
0.186 


0.527 
0.137 


0.529 
0.134 


0.620 
0.135 



We would like to note that for the modified confining potential the constituent mass 
grows by only about 15% for the D-wave states and about 25% for the P-wave states, in 
contrast to the situation for the standard linear potential (see Table |TT|) where this growth 
is substantially larger. Also the matrix elements < > for the nD states with rir = 1,2,3 
turn out to be equal within 2%. 
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